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ABSTRACT 

In this paper, we report the detection of ejecta in the middle-aged Galactic super- 
nova remnant G296. 1—0.5 with the X-ray Imaging Spectrometer onboard the Suzaku 
satellite. The spectra of three lobes, north, southeast and southwest and inter-lobe 
regions, consist of soft (0.3—2.0 keV) emission originated from non-equilibrium ion- 
ization plasma. In north, southeast and inter-lobe regions, the thermal emission can 
be represented by a one-component, in southwest region it can be represented by two- 
component non-equilibrium ionization (VNEI) model. The spectra of studied regions 
have lines of N, O, Ne, Mg and Si elements. Si emission from this remnant is shown for 
the first time in this work. Enhanced abundances of Ne, Mg and Si elements obtained 
show the ejecta contribution in all regions. Assuming that the remnant is in Sedov 
phase, we obtained ambient density hq ~ 0.45 cm~ 3 , age t ~ 2.8 x 10 4 yr, shock 
velocity V s ~ 320 km s _1 , shock temperature T s ~ 1.2 x 10 6 K, and swept-up mass 
M sw ~ 340 M© at an adopted distance of d=3 kpc. 
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1 INTRODUCTION 

The shell type Galactic supernova remnant (SNR) 
G296.1-0.5 (RA(2000) = ll h 51 m 10 s , Dec. (2000) = 
— 62°34') with 37x25 arcmin 2 angular size was discovered 
in rad i o band at 408 MHz (Clark, Caswell fc Green ||1973 



|1975[ ). |Longmore, Clark fc Murdin| ( |1977[ ) performed opti- 
cal observation of the remnant and reported the identifi- 



cation of the nebulosity and its filamentary nature. |Hutch-| 
|ings, Crampton fc Cowley] ( |1981[ ) suggested that there was 
another SNR located at the northeast part of the radio rem- 
nant. Since the region in which the remnant is located on 
the Galactic plane is rich in SNRs, this suggestion seemed 



reasonable. Caswell & Bornes (19831, with an improved ra- 



dio map at 408 MHz, and Bignami et al. (19861, with X 



ray data from EXOSAT and Einstein observations, reported 
that this is a single, complex remnant. From ROSAT X-ray 



data Hwang & Market ( 1994 I obtained an interstellar ab- 
sorbtion of Ah ^ 1.5 x 10 cm -2 and a low electron tem- 
perature kT c ~ 0.2 keV from single temperature model fit- 
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ting. They also applied two temperature component model, 
obtained electron temperatures of ~0.1 keV and ~3.5 keV 
for each component. Assuming the distance of 4 kpc and el- 
lipsoidal shell geometry for the remnant, they found total X- 
ray emitting mass and ambient density to be ~250 M0 and 
~0.8 cm -3 , respectively. Using the Sedov equations, they 
found an age of ~2 x 10 4 yr. From the data of XMM-Newton 



observation, Castro et al. ( 2011 1 found the electron temper- 
ature fcT e =0.6 keV by applying a non-equilibrium ionization 
model (VNEI) to the spectra. They obtained a low absorb- 
ing column density, Ah = (2 — 4) x 10 20 cm -2 , excess of 
N, and a deficit of O. They suggested a massive progenitor 
star and concluded that the remnant was resulted from core- 
collapse supernova. They have detected a transient source 
2XMMi Jl 15004.8622442 at the edge of this remnant. 

Various distance estimates were given from different 



studies; Longmore, Clark & Murdin ( 1977 1 estimated the 
distance d=3 ± 1 kpc by using Ha? and H/5 observations 



and reddening measurements in the direction of crux ( |Miller 
19721, and by fitting the radial velocity (—35 km s~ with 
respect to the local standard of rest) to galactic rotation. 



From E— D relation Caswell & Bornes ( 1983 1, Case & Bhat 



tacharya ( 1998 1 and Clark, Caswell & Green 



(1973) adopted 



d=7.7 kpc, d=6.6 kpc and d=4.9 kpc, respectively. |Castro| 



et al. (20111 assumed d=2 kpc, considering the giant HI 
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shell, GSH 304-00-0.5 located in front of the remnant at d 
~ 1.2 kpc. 

Our aim is to study middle-aged SNRs to investigate 
interstellar medium (ISM) /shock interactions and ISM it- 
self, i.e., the status of the ionization state, temperature (line 
emission) variation all over the remnant. G296. 1—0.5 is par- 
ticularly interesting one with its complex shock and plasma 



structure. The Japanese X-ray observatory Suzaku (Mitsuda 
et al.|2007[ |, since it has a large collecting area and low back- 
ground, is capable of resolving line emission characteristic 
of SNRs, especially at low energies, that provides valuable 
information in understanding SNRs. By using the archival 
data of Suzaku, we were able to produce higher quality im- 
age and the spectra of the remnant, which motivated us for 
this study. 

This paper is organized as follows. We describe the ob- 
servations and the data reduction in Section 2. In Section 
3, we explain the details of image and spectral analysis. In 
Section 4, we discuss our results and description of the X-ray 
emission from G296. 1—0.5. 



2 OBSERVATIONS AND DATA REDUCTION 

The Suzaku observations of the north and south regions of 
G296.1-0.5 were made on 2007 August 09 and on 2008 
January 17, respectively, with X-ray Imaging Spectrometer 



(XIS: Koyama et al. ( 2007 1). The observation ID and expo- 
sure time are 502068010, 502069010 and 77.2, 69.2 ksec, re- 
spectively. The XIS consists of four sets of X-ray CCD cam- 
era system (XIS0, 1, 2, and 3). XIS1 has a back-illuminated 
(BI) sensor, while XIS0, 2, and 3 have front-illuminated 
(FI) sensors. The XIS2 sensor was available only until 2006, 
therefore, we use data of XIS0, XIS1, XIS3. The XIS has a 
field of view (FOV) of a 17.8x17.8 arcmin 2 . 

In both observations, the XIS was operated in the nor- 
mal full-frame clocking mode with the standard 3x3 and 
5x5 editing mode. Response matrices and ancillary re- 
sponse files (ARFs) were generated for each XIS indepen- 
dently using xissimrmfgen version 2007-05-14 and XISSI- 
MARFGEN version 2008-04-05 ( jlshisaki et al.||2007| . Reduc- 
tion and analysis of the data were performed following the 
standard procedure using the HEADAS v6.4 software pack- 
age, and spectral fitting was performed with xspec v. 11.3.2 
( Arnaud|1996] ). 



3 ANALYSIS 

3.1 Image Analysis 

Figure la and lb present the smoothed XIS0 images of north 
and south regions of G296. 1—0.5 in 0.3—10 keV energy band, 
respectively. In the northern part there is one bright lobe, 
while in the southern part there are two bright lobes. So, to 
determine the temperature variations and ionization states, 
we studied the bright lobes and a region between the lobes 
in the south part as shown in Fig. la and lb with solid el- 
lipses. The regions are abbreviated as north (N), southeast 
(SE), southwest (SW) and inter-lobe (M) region. M region 
is selected as the reference location in order to highlight 
the peculiar properties of the lobes if there is any. Dashed 
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Figure 1. Suzaku XIS0 image of (a) N region, (b) SE, SW and 
M regions of G296. 1—0.5 in the 0.3—10 keV energy band (gray 
scale). Regions for the source and the background are indicated 
with the solid and dashed ellipses, respectively. Two corners of 
each FOV, where the calibration source of 55 Fe is illuminated, 
are masked in these images. The coordinates (RA and Dec.) are 
referred to epoch J2000. 



ellipses in both figures show the regions chosen for back- 
ground subtraction. The calibration sources located at the 
two corners in the FOV are excluded. 

Figure 2 gives the mosaic image of XIS0 in 0.3—10 keV 
energy band, N and S, which are overlaid with the radio 
image obtained at 843 MHz by |Whiteoak fc Green] (|1996[) 
for comparison, the positions of optical filaments ( Longmore, 
Clark fc Murdin|1977[ ), Ha emission and thin filament (|Rus- 
seil & Parker 2001) and molecular material ( |Brand et al. 
19871 are also pointed out. 
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Figure 2. Suzaku X-ray mosaic image (blue for N and green 
for SE and SW) and its radio image (red) obtained at 843 MHz 
are overlaid and smoothed with a Gaussian kernel of <r=2 arcsec, 
and units are in counts s — 1 . The linearly spaced radio contours 
are 0.02, 0.06, 0.10, 0.14, 0.18 and 0.22 counts/pixel. The blue 
and green squares indicate the XIS FOV for N and S regions, 
respectively. The positions of optical filaments ||Longmore, Clark 
fc Murdin|[l977[ l, Ho emissio n, thin filament l|Russeil fc Parker 
2001[ | and molecular material | |Brand et al.|1987[ l are marked as, 
plus, circle, cross and diamond, respectively. The coordinates (RA 
and Dec.) are referred to epoch J2000. 



3.2 Spectral Analysis 

We extracted the representative XIS spectra from four ellip- 
tical regions named as N, SE, SW and M (Fig. la and lb) by 
using xselect version 2.4. The coordinates of the ellipses, 
their sizes and angles are listed in Table 1. All spectra are 
grouped with a minimum 20 counts bin -1 . 

From the spectra we see that emission in 0.3—2.0 keV 
soft energy range is dominant and K-shell lines of N, O, 
Ne, Mg, and Si are present for all regions, while in 2—10 
keV energy range the emission is too faint. So, we select the 
energy range of 0.3—2.0 keV for spectral fitting. 

To analyze XIS data of the selected regions, since previ- 
ous XMM-Newton observations showed non-equilibrium ion- 
ization (NEI) condition in this SNR, we applied NEI colli- 
sional plasma with variable abundances model (VNEI model 
in XSPEC, Borkowski, Lyerly & Reynolds ( 2001 1) modified 
by interstellar absorption (wabs model in XSPEC, |Morri- 



son & McCammon ( 1983 1) for all spectra. For N, SE and 
M regions, one-component VNEI model best describes the 
spectral characteristics of the emission. In the SW region, 
a two-component VNEI model fits (xV d -°- f -=349. 6/321) 
the emission better than a one-component VNEI model 
(x 2 /d.o.f.=403. 6/324) statistically. During the analysis, to 
see if there are small scale variations the absorbing col- 
umn density (-/Vh), electron temperature (kT a ), the ioniza- 
tion timescale (r, which is often used as a key diagnostic of 
the NEI state, is the product of the electron density mul- 



tiplied by the time after the shock heating) are allowed to 
vary. Since the emission lines of N, O, Ne, Mg, and Si are 
clearly seen in the spectra, their abundances are set as free 
parameters while the other elemental abundances are fixed 



at their solar values (Anders & Grevesse 19891. The best- 



fitting parameters, their errors and % values for all regions 
are given in Table 2 with corresponding errors at the 90 per 
cent confidence limit. We note that, due to better photon 
statistics of BI (XIS1) than FI (XIS0 and XIS3), only XIS1 
results are given in this table. The background-subtracted 
XIS1 spectra in the 0.3-2.0 keV energy band for N, SE, SW 
and M regions are shown in Figure 3. We also plot the confi- 
dence contours of the absorbing column density (Nh) versus 
the electron temperature (kT e ) and N abundance for all re- 
gions in Figure 4. Both parameters are very well constrained 
for N and SE regions. 



4 DISCUSSION AND CONCLUSIONS 

In this paper, we give a description of the X-ray emis- 
sion from G296. 1—0.5 based on Suzaku archival data. We 
obtained a clear image and high quality spectra of X-ray 
emission that lead us detection of ejecta component in this 
middle-aged remnant. 

The combined image, Fig. 2, shows that the remnant ex- 
tends to a large area (37x25 arcmin 2 ) in radio band, while 
it is concentrated in three regions as bright lobes in X-ray 
band, two being in southern part and one being in the north- 
ern part. Suzaku observed the lobes shown by squares, not 
the entire remnant. Therefore we could study only the lobes 
and a region between SE and SW lobes. The bright X-ray 
emission from lobes is well correlated with radio emission. 
In both bands, G296. 1—0.5 shows partial (irregular) shell 
morphology. Since radio emission is not effected by absorb- 
tion we can say that the shape of the remnant is intrinsically 
anisotropic and it is interacting with a non-uniform medium. 
The optical imaging performed byjLongmorc, Clark & Mur- 
|dhT| ( fl977l ) shows that SW region covers optical filaments 
with [S ll]/H« ~ 1.2 which is a typical value for shocked gas. 
They obtained an electron density as high as 2 x 10 3 cm -3 . 



Russeil & Parker (20011 detected strong localized diffuse 



Ha emission using the data from the Marseille survey and 
many thin filaments by equivalent AAO/UKST Ha survey 



in the direction of the remnant. Brand et al.| ( |1987| ) detected 
molecular material with a radial velocity of ~ —38 km s _1 
at RA(2000) = ll h 50 m 27 s , Dec. (2000) = -62°38'55" by 
using CO observations (see Fig. 2). All these observations 
may indicate that a shock-cloud interaction is going on and 
the lobes may result from these interactions. 

X-ray emission of regions N, SE and M, can be described 
by one-component non-equilibrium ionization model (VNEI) 
in 0.3—2.0 keV band with electron temperature kT e in the 
range of (0.51 — 0.76) keV which is in good agreement with 
the recent work by |Castro et al. ](|2011[), wh ile almost three 
times higher than that of Hwang & Market ( 1994 1. The ion- 
r, ranges between (1.25 



ization time scale, r, ranges between (1.25 — 2.46) x 10 10 
~ 3 s indicating that the plasma in these regions is far 
from ionization equilibrium condition (r is typically re quired 
to be JSlO 12 cm~ 3 s for full ionization equilibrium (Masai 
1 1984[ >) . On the other side, in the SW region, thermal emis- 



sion can be adequately described by two-component non- 
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equilibrium ionization model (VNEI+VNEI): the colder, 
fcT c ~ 0.18 keV, is characterized by lower solar abundances 
of O and N and ionization time scale r ~ 5.4 x 10 11 cm _3 s, 
the hotter, kT c ~ 0.84 keV, by over abundances of Ne, Mg 
and Si and ionization time scale r ~ 1.7 x 10 10 cm _3 s, both 
components are far from ionization equilibrium condition in 
this region also. 

The elemental abundance values of N, O, Ne, Mg and Si 
show small variations in N, SE and M regions. In N and SE 
regions, the abundance of O appears to be lower or about so- 
lar value, N, Ne, Mg and Si values are somewhat higher than 
solar. In general, shocked ISM spectra shows lower or about 
solar abundances. The plasma in these two regions contains 
swept-up ISM and ejecta contamination. In M region, abun- 
dance values and their errors seem somewhat higher as com- 
pared to N and SE regions (see Table 2). We note that the 
photon statistics in this region is poor (see Fig. 3), there- 
fore we can say that, like N and SE regions, the plasma 
contains swept-up ISM and ejecta contamination. In SW re- 
gion, the abundances of N and O from cold component are 
lower solar indicating ISM origin, while Ne, Mg and Si values 
from hot component are high enough to show the evidence of 
ejecta origin in this region also. G296. 1—0.5 is a middle-aged 
remnant, thus the presence of ejecta at its evolved stage is 
not expected. However, recent X-ray observations show that 
in a number of middle-aged remnants metal rich ejecta are 
still present (e.g. G156.2+ 5.7 (|Yamauchi et al.|1999[ ), N49B 
( Park et al.||2003[), IC443 (|Troja et al.||2008| > and Puppis A 
(Hwang, Petre & Flanagan||2008l)). As seen in Fig. 2, the 



presence of molecular material in the vicinity of SE and SW 
regions may cause the survival of ejecta in especially these 
regions of this middle-aged remnant. 

The absorbing column density Nn value is ranging 
between (2.0 — 14.6) xlO 20 cm -2 across the remnant. Ah 
value of SW region obtained from two temperature model, 
~14.6 x 10 20 cm -2 , is in agreement with the A^h value of 



Hwang & Market (1994) obtained by applying two temper- 



ature model. On the other hand, A^h value of N and SE 
regions are 2 — 3 times smaller than A^h value of |Hwang fc| 
Market ( 1994[ ) obtained by applying one temperature model, 



while 4 — 5 times higher than that of Castro et al. (2011 1 



Therefore, for our calculations we will use d=3 kpc value, 
that is in between the values adopted by |Hwang fc Market| 
( 1994 1 and Castro et al. ( 2011 1 and is in agreement with the 
value estimated by |Longmore, Clark fc Mu7dTn| ( [19771 ). We 
note that best-fitting A^h value for SW region is relatively 
higher than the other regions (see Table 2). The reason for 
this high Ah value may be the inhomogeneity of the fore- 
ground absorbing medium (see Fig. 2). Since an accurate 
estimate of Ah put constraints on the distance to the rem- 
nant, we obtained the confidence contours of Nu versus kT c 
and N abundance as given in Fig. 4. We see that for N and 
SE regions both parameters are very well constrained. The 
X-ray flux and the corresponding luminosity values are ob- 
tained to be F x ~ 1.1 xlO" 9 ergs _1 cm" 2 and L x ~ 1.2xl0 36 
erg s _1 at d=3 kpc in the 0.3—2.0 keV energy range. 

It is suggested that most optically observed SNRs are 



in the adiabatic phase of their evolution (McKee & Cowie 



1975). Furthermore the large physical size of this remnant, 
about 46 pc in radio, also suggests that it is in adia- 
batic expansion phase. We calculated physical parameters 



of G296. 1—0.5 by using information obtained by XIS im- 
ages from Sedov equations ( Sedov| 19*59 1 : 

t ,2/5, E 



R 



4 x 10 19 ( 



10 4 yr 



) 2/5 ( 



) 1/5 (: 



10 51 erg 1cm 
E 



) 1/5 (: 



no 



-1/5, 



10 51 erg 1cm 



(1) 



-\(2) 



E ^2/5/ n N-2/S-, 



T - 3Xl06 (l(fe)" 6/5 (ll^)^&)"^ K ' ^ 



L x = 3.3 x 10 29 n 2 Tc°- 5 i? s 3 ergs" 1 
4 3 

Mism = -nR s (imnno, 



(4) 
(5) 



where Ra is the shock radius, ~23 pc, the explosion en- 
ergy is assumed to be Esi = l, \i is atomic weight of 0.604, toh 
is mass of a hydrogen atom. Through the equations (1—5), 
we obtained ambient density no ~ 0.45 cm -3 , age of the 
remnant, t ~ 2.8 x 10 4 yr, confirming that G296. 1—0.5 is 
a middle-aged SNR. The corresponding ionization age n e t 
is then ~2.7 x 10 11 cm~ 3 s, a plasma with this ionization 
age is in non-equilibrium ionization condition. Assuming 
the electron-ion temperature equilibration, shock velocity 
V s that corresponds to best-fitting electron temperature be- 
hind the main shock front (kT c ~ 0.18 keV) is calculated to 
be ~320 km s _1 . Assuming uniform density for the ambi- 
ent ISM, we estimated shock temperature T B ~ 1.2 x 10 6 K 
and the mass swept-up by the blast wave M aw ~ 340 M0. 
Castro et al. (20111 suggested that G296. 1—0.5 is resulted 
from a core-collapse supernova with a massive (> 25 M0) 
progenitor star. We compared our best-fitting abundances 
of N, O, Ne and Mg relative to Si for SE, SW and M regions 
with abundance ratios in core-collapse model of Wooslcy & 
Weaver (1995]) and the obtained results are shown in Figure 
5. As seen from the figure, we can say that the progenitor 
star mass might be > 25 M©. Considering the high N abun- 
dance value it appears to be roughly consistent with 30 MQ. 
When we compare this value with our estimated swept-up 
mass which is approximately ten times higher, we can say 
that G296. 1—0.5 is at its critical stage and the contamina- 
tion of ejecta in this middle-aged SNR is not surprising. 
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Table 1. Centre coordinates, sizes and angles of N, SE, SW, M and background regions. 



Regions 


Centre Coordinates 


Size 


Angle 




RA(2000) ; Dec. (2000) 


(arcmin) 


(degrees) 




^ h m s 


. o / //^ 






N 


11 50 59 


-62 19 46 


7.5 X 3.0 


318 


Background" 


11 51 45 


-62 25 56 


2.4 X 1.8 





SE 


11 51 12 


-62 39 13 


2.7 x 6.3 


332 


SW 


11 50 02 


-62 44 38 


2.8 x 5.2 


339 


M 


11 50 20 


-62 37 55 


1.9 x 3.5 


334 


Backgrounoj^] 


11 51 37 


-62 48 19 


1.4 x 1.9 






a for N region; 
^ for S region. 



Table 2. Best spectral fitting parameters of XIS1 spectra obtained for N, SE, SW and M region of G296. 1—0.5 in 0.3- 
corresponding errors at 90 per cent confidence level. 



2.0 keV with 



Component 


Parameters 


N 


SE 


SW 


M 


Absorbtion 


A^xlO 20 cm~ 2 ) 


7.5 ±0.5 


4.7 ±0.6 


14.6 ±1.4 


2.0 ±0.5 


VNEI 


AiT c (keV) 


0.58 ±0.02 


0.51±0.02 


0.18±0.01 


0.76±0.15 


Abundance]^] 


N 


1.91 ±0.07 


2.15 ±0.09 


0.75 ±0.12 


3.08 ±0.36 





0.81 ±0.02 


1.01 ±0.03 


0.39 ±0.03 


0.91 ±0.07 




Ne 


1.56 ±0.04 


1.92 ±0.05 


(1) 


2.3±0.2 




Mg 


1.37±0.06 


1.44±0.10 


(1) 


2.3±0.5 




Si 


1.26±0.21 


1.75±0.36 


(1) 


2.7±1.4 




n c i(xl0 10 cm" 3 s) 


2.24 ±0.15 


2.46±0.21 


54.2±1.3 


1.25±0.29 




e.mE 


5.8±1.2 


4.6±1.1 


89.3±10.2 


0.59±0.03 


VNEI 


A:T (keV) 






0.84±0.09 






Ne 






2.38±0.31 






Mg 






1.55±0.21 






Si 






2.39±0.55 






n c i(xl0 10 cm" 3 s) 






1.68±0.14 






E.M. 






1.8±0.4 






X 2 /d.o.f. 


619.2/380=1.63 


590.2/338=1.75 


349.6/321=1.09 


202.7/185=1.09 



Abundance ratio relative to the solar value I Anders fc Grevesse|198"9 i, (1) indicates that the elemental abundance is fixed at solar 



Emission measure EM= J n e nndV in the unit of 10 s ' cm a , where n e and njj are number densities of electrons and protons, respectively 



and V is the X-ray-emitting volume. 
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Figure 3. XIS1 spectra of G296. 1—0.5 obtained in the 0.3—2.0 keV energy band. The indicated fits consists of an absorbed VNEI 
model for N, SE and M regions and an absorbed VNEI+VNEI model for SW region. The bottom window gives the residuals from the 
best-fitting model for XIS1 spectra. 
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Figure 4. Confidence contours of TVh versus kT e (left) and Nh versus N abundance (right) spectral fitting for all regions of G296. 1—0.5. 
The contours are at the 68, 90 and 99 per cent confidence levels. The parameter values corresponding to the best-fitting are marked by 
a plus sign (+). 
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Figure 5. Best-fitting abundances of N, O, Ne and Mg relative to Si relative to solar (for SE, SW and M regions) are compared with 
the theoretical predictions of core-collapse model ||Woosley &; Weaver|1995[l with 25, 30 and 35 M0 progenitors. 



